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Abstract 

The  electrochemical  behavior  of  a  direct  ethanol  feed  proton  exchange  membrane  fuel  cell  (DEFC)  operating  under  steady-state  isothermal 
conditions  at  1  atm  at  both  anode  and  cathode  sides  is  considered.  A  mathematical  model  that  describes  in  one  phase  and  one  dimension  the  ethanol 
mass  transport  throughout  the  anode  compartment  and  proton  exchange  membrane  is  developed.  The  influence  of  the  operation  parameters  such 
as  current  density,  temperature,  catalyst  layer  thickness  and  ethanol  feed  concentration  on  both  anode  overpotential  and  ethanol  crossover  rate  has 
been  examined.  According  to  the  simulation  results,  it  was  found  that  the  anode  overpotential  is  more  sensitive  to  the  protonic  conductivity  than 
to  the  diffusion  coefficient  of  ethanol  in  the  catalyst  layer.  It  was  concluded  that  in  the  case  of  low  current  density  values  and  high  concentrations 
of  ethanol  aqueous  solutions,  ethanol  crossover  is  a  serious  problem  for  a  DEFC  performance.  Finally,  it  was  found  a  good  agreement  between 
simulation  and  experimental  results. 

©  2006  Elsevier  B.V.  All  rights  reserved. 

Keywords:  Mathematical  modeling;  Direct  ethanol  fuel  cell  (DEFC);  Ethanol  crossover 


1.  Introduction 

The  last  decade  direct  ethanol  fuel  cells  (DEFCs)  have 
attracted  more  and  more  interest  due  to  ethanol  intrinsic  advan¬ 
tages  [1-21].  It  has  already  been  recognized  that  ethanol  is 
a  promising  fuel  which  is  non-toxic,  renewable,  and  can  be 
easily  produced  in  great  quantity  by  the  fermentation  of  sugar- 
containing  raw  materials  [1,2].  Moreover,  a  high  theoretical 
mass  energy  density  (about  8.00  kWh  kg~ 1 )  [  16]  provides  it  with 
a  potential  candidate  fuel  for  proton  exchange  membrane  fuel 
cells.  However,  there  are  some  challenges  that  DEFCs  meet,  such 
as  the  high  anode  overpotential  values,  slow  ethanol  electro¬ 
oxidation  kinetic  and  fuel  crossover  from  the  anode  to  cathode 
through  the  membrane,  which  inevitably  decrease  the  fuel  cell 
efficiency.  At  the  same  time,  the  permeated  ethanol  and  its  oxida¬ 
tion  intermediate  products  could  poison  the  cathode  catalyst.  In 
our  previous  work  the  experimental  results  of  ethanol  crossover 
have  already  been  reported  [9]. 

Mathematical  modeling  is  essential  for  the  development  of 
fuel  cells  because  it  allows  extensive  comprehension  of  the 
parameters  affecting  the  performance  of  a  single  fuel  cell  or 
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fuel  cell  systems  [23-26].  The  work  of  Haraldsson  and  Wipke 
[3]  is  directed  to  the  analysis  and  evaluation  of  fuel  cell  system 
models.  Jeng  and  Chen  [12]  developed  a  mathematical  model  in 
order  to  predict  the  behavior  of  a  direct  methanol  fuel  cell  and 
examine  the  effect  of  the  current  density  on  methanol  crossover. 
However,  there  are  some  empirical  models  [8]  developed  specif¬ 
ically  for  various  applications  and  operating  conditions. 

In  the  present  work,  a  mathematical  model  for  the  simulation 
of  a  direct  ethanol  fuel  cell  anode  has  been  developed  with  the 
purpose  to  further  investigate  the  effect  of  fuel  cell  main  param¬ 
eters  (current  density,  temperature,  ethanol  concentration)  on 
both  ethanol  crossover  and  anode  overpotential.  More  precisely, 
this  model  aims  at  the  solution  of  the  governing  equations  of 
mass  transport  through  both  the  catalyst  layer  and  the  proton 
exchange  membrane. 

2.  Theory 

As  schematically  shown  in  Fig.  1,  the  ethanol  fed  to  the 
anode  compartment  of  a  direct  ethanol  fuel  cell  could  be  par¬ 
tially  reacted,  partially  unreacted  and  as  the  experimental  results 
[9,10]  showed,  partially  crossovered  to  the  cathode  compart¬ 
ment.  It  can  also  be  seen  that  the  anode  compartment  consists  of 
diffusion  layer  and  catalyst  layer.  The  diffusion  layer  is  made  of 
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Nomenclature 

A„ 

electrode  area  where  the  reaction  takes  place 

CEtOH 

local  ethanol  concentration 

/^>an 

^EtOH 

ethanol  concentration  on  the  anode  side  (M) 

/'-'cath 

^EtOH 

ethanol  concentration  on  the  cathode  side  (M) 

CpEtOH 

feed  ethanol  concentration  (M) 

/^ref 

^EtOH 

reference  ethanol  concentration  (M) 

Cs,EtOH 

ethanol  concentration  at  the  surface  of  diffusion 
layer  (M) 

r° 

'-EtOH 

ethanol  concentration  at  z  =  0  (M) 

nd,eff 

^EtOH 

effective  diffusion  coefficient  at  the  diffusion 
layer  (cm2  s-1) 

nc,eff 

^EtOH 

effective  diffusion  coefficient  at  the  catalyst  layer 
(cm2  s-1) 

T^d.EtOH 

diffusion  coefficient  of  ethanol  in  water  (cm2  s-  1 ) 

nm 

^EtOH 

diffusion  coefficient  of  ethanol  in  PEM  (cm2  s-  1 ) 

nm 

uh2o 

diffusion  coefficient  of  water  in  PEM  (cm2  s-  1 ) 

F 

Faraday  constant  (96498  C  mol-1) 

i 

protonic  current  density  (A  cm-2) 

*o,ref 

reference  exchange  current  density  (A  cm-2) 

i 

cell  current  density  (A  cm-2) 

mass  transfer  coefficient  of  ethanol  from  the  feed 
stream  to  diffusion  layer  (cm  s-1) 

kd 

mass  transfer  coefficient  of  ethanol  in  diffusion 
layer  (cm  s-1) 

km 

mass  transfer  coefficient  of  ethanol  in  PEM 
(cms-1) 

Km 

protonic  conductivity  of  ionomer 

Jy'eft 

Km 

effective  protonic  conductivity  in  catalyst  layer 

Kf 

effective  conductivity  of  solid  phase  in  catalyst 
layer 

KS 

electronic  conductivity  of  solid  phase  (Pt-Ru/C) 

(S  cm-1) 

k 

thickness  of  catalyst  layer  (cm) 

k 

thickness  of  diffusion  layer  (cm) 

PEM  thickness  (cm) 

Mn2o 

molecular  weight  of  water  (g  mol- 1 ) 

AfetOH 

local  ethanol  flux  in  catalyst  layer  (mol  cm-2  s- 1 ) 

-,vEtOH 

ethanol  flux  through  the  diffusion  layer 
(mol  cm-2  s-1) 

Nm 

JVEtOH 

ethanol  flux  through  PEM  (mol  cm2  s- 1 ) 

iVdiff 

water  flux  due  to  the  diffusion  mechanism  in  PEM 
(mol  cm-2  s-1) 

A^eiectr  drag  water  flux  due  to  the  electro-osmotic  drag  in 
PEM  (mol  cm-2  s-1) 

Ah2o 

local  water  flux  in  catalyst  layer  (mol  cm-2  s- 1 ) 

Nd,n 2o 

water  flux  through  the  diffusion  layer 
(mol  cm-2  s-1) 

Nm 

-,vh2o 

water  flux  through  PEM  (mol  cm-2  s-1) 

At 

total  flux  of  water  and  ethanol  in  catalyst  layer 
(mol  cm-2  s-1) 

R 

universal  gas  constant  (8.3144  J  mol-1  K-1) 

T 

cell  temperature  (K) 

ud 

superficial  velocity  of  water  in  diffusion  layer 
(cms-1) 

vm 

superficial  velocity  of  water  in  PEM  (cms-1) 

jEtOH 

mole  fraction  of  ethanol 

Greek 

letters 

«a 

anode  transfer  coefficient 

£d 

void  fraction  of  diffusion  layer 

£C 

void  fraction  of  catalyst  layer 

Y 

order  of  reaction 

il 

anode  overpotential  (V) 

(Pm 

potential  of  the  ionomer  phase  (V) 

potential  of  electronic  conduction  phase  (V) 

electro-osmotic  drag  coefficient  of  water 

Ph2o 

density  of  water  (1.0  gem-3) 

a  porous  and  electrically  conductive  material,  through  which  the 
electrons  generated  in  the  anode  catalyst  layer  are  transported  to 
the  current  collector.  On  the  other  hand,  the  catalyst  layer  is  the 
place  where  the  reaction  of  the  ethanol  electro-oxidation  takes 
place,  releasing  protons  and  electrons.  Both  protons  and  elec¬ 
trons  are  then  transported  to  the  cathode:  the  electrons  through 
an  external  circuit  while  the  protons  through  a  proton  exchange 
membrane  (PEM). 

The  mathematical  model  has  been  developed  taking  into  con¬ 
sideration  the  following  assumptions:  (a)  the  cell  operates  under 
steady-state  isothermal  conditions,  (b)  the  pressure  in  both  anode 
and  cathode  compartment  is  1  atm,  (c)  the  model  equations  are 
defined  in  one  dimesion  and  (d)  the  origin  is  set  at  the  interface 
between  the  diffusion  layer  and  the  catalyst  layer. 

2.1.  Flow  channel 

The  ethanol  concentration,  which  is  fed  into  the  flow  chan¬ 
nel,  is  indicated  as  feed  concentration  (Cf,ei:Oh)-  The  ethanol 
transportation  from  the  feed  stream  to  the  diffusion  layer  can  be 
described  by  the  following  equation: 

^EtOH  =  £g(CF,EtOH  —  C^EtOH)  (1) 

where  AgtOH  stands  for  ethanol  flux  through  diffusion  layer, 
Cp.EtOH  ethanol  feed  concentration,  CSjEtOH  ethanol  concentra¬ 
tion  at  the  surface  of  the  diffusion  layer,  and  kg  mass  transfer 
coefficient.  The  mass  transfer  coefficient  from  the  feed  stream 
to  the  diffusion  layer  can  be  approximately  expressed  by  using 
the  empirical  equation  [12]: 

/  /  \ 032 

k„  =  1.87  x  10-4 -  (2) 

s  \  0.003  ) 

where  I  denotes  the  current  density  of  the  fuel  cell. 

2.2.  Diffusion  layer 

Diffusion  layer  is  made  of  carbon  particles  that  are  inactive  to 
the  electrochemical  reaction  for  a  DEFC.  After  passing  through 
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Fig.  1.  Schematic  presentation  of  a  DEFC  and  ethanol  transportation  through  anode  compartment  and  PEM. 


diffusion  layer,  water  is  partially  consumed  in  the  electrochemi¬ 
cal  reaction  which  takes  place  in  the  catalyst  layer.  The  remain¬ 
ing  quantity  could  further  migrate  through  PEM  by  diffusion, 
electro-osmosis  and  hydraulic  permeation.  The  electro-osmotic 
drag  appears  when  the  electrochemical  cell  operates,  while  the 
hydraulic  permeation  appears  in  the  presence  of  pressure  dif¬ 
ference  between  anode  and  cathode.  As  above  mentioned  the 
pressure  has  been  assumed  at  both  anode  and  cathode  sides  equal 
to  1  atm  and  consequently  the  water  flux  through  the  diffusion 
layer  can  be  attributed  to  both  diffusion  and  electro-osmosis 
phenomena  and  described  by  the  following  equation: 

^d.H20  =  ^  +  ^H20  (3) 

where  F  is  the  Faraday  constant  and  /V|"[,i()  the  water  flux  through 
PEM. 

Through  the  diffusion  layer  ethanol  is  transported  by  diffu¬ 
sion  and  convection: 

j»d  nd,eff  dCaOH  , , 

^EtOH  =  -  ^EtOH  — - h  VC.EtOH  (4) 

where  is  the  effective  diffusion  coefficient  of  ethanol  in 

the  diffusion  layer,  vc.EtOH  the  local  mole  fraction  of  ethanol  in 
the  diffusion  layer,  and  Nj  the  total  flux  of  ethanol  and  water. 
The  effective  diffusion  coefficient  through  a  porous  media,  the 
bulk  diffusion  coefficient  (/Jd.EtOii)  and  the  void  fraction  of  the 
diffusion  layer  (ed)  are  related  by  Bruggeman’s  empirical  equa¬ 
tion  [19,20], 


The  local  mole  fraction  of  ethanol  yc.htoi  I  is  given  as  follows: 

ChtOI  I  CEtOH  /Wh2oChk)|[ 

yc.EtoH  =  - - — - =  — - = - - -  (o) 

<-H20  +  L-KlOH  Ch20  PH20 

Combining  Eqs.  (6)  and  (4),  Eq.  (4)  can  be  expressed  by  a 
first  order  differential  equation  as  follows: 


-^EtOH  —  ~D 


d,eff  dCEtOH  , 
Eton  d,  + 


M^oCeiow  .. 

- - - A^d.H2o 

Ph2o 


(7) 


Assuming  that  /^Eton  and  ph2o  are  constant,  Eq.  (7)  can  be 
solved  within  the  interval  of  —  ld<z<0  taking  the  following 
form: 


-^EtOH  — 


Cs,EtOHe1’  /k  —  C 

^/FT T“ 


0 

EtOH  L,d 


(8) 


where  kd  =  D^f^/ld  is  the  mass  transfer  coefficient  in  the  dif¬ 
fusion  layer,  CgtOH  the  ethanol  concentration  at  the  interface 
between  the  diffusion  layer  and  the  catalyst  layer  (i.e.  at  -  =  0). 
i>d  =  Mh2o  A^h2o/Ph2o  denotes  the  superficial  velocity  of  water 
in  the  diffusion  layer  and  Id  the  thickness  of  the  diffusion  layer. 

By  combining  Eqs.  (1)  and  (8)  and  eliminating  Cs>EtOH,  the 
expression  of  the  ethanol  flux  through  the  diffusion  layer  as 
a  function  of  the  ethanol  feed  concentration  (Cf.eioh)  can  be 
obtained: 


-^EtOH  — 


Cf.EtOHe1  ,k  -  CgtOH 
svi/ki[(vd/kd)+  1]  -  1 


(9) 


2.3.  Catalyst  layer 


D 


d.eff 

EtOH 


-fpdt3/2n 

—  (£  >  p^d.EtOH 


Catalyst  layer  is  the  place  where  the  reaction  of  ethanol 
(5)  electro-oxidation  takes  place.  The  rate  of  the  electrochemical 
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reaction  can  be  described  by  using  the  Butler- Volmer  equation 
and  in  a  simpler  way  by  the  Tafel  equation: 


di 

dz 


^ido,ref 


CkiOI  I  \ 

rref  / 

'-'EtOH  / 


Y 

exp 


^12«a Frf^ 


(10) 


where  i  is  the  local  protonic  current  density,  Av  the  specific 
area  of  the  reaction  surface,  io  ref  the  reference  exchange  cur¬ 
rent  density,  y  the  order  of  reaction,  CEtOH  the  local  ethanol 
concentration  in  the  catalyst  layer,  Cj^fon  the  reference  ethanol 
concentration,  which  is  associated  with  f0,ref ,  a a  the  anode  charge 
transfer  coefficient  and  rj  is  the  anode  overpotential. 

The  driving  force  for  both  water  and  ethanol  transporta¬ 
tion  through  the  catalyst  layer  consists  of  the  diffusion  and  the 
electro-osmotic  drag. 


iVEtOH  —  —D 


c.eff  dCEtOH 
EtQH  dz 


A^H2oCEtOH  .. 

- - - Ah20 

Ph2o 


(ID 


c  eff 

where  (  is  the  effective  diffusion  coefficient  of  ethanol  in 
the  catalyst  layer,  and  Ah20  is  the  local  water  flux.  The  water 
mass  balance  gives: 


Nn2o  —  +  nv\20  (12) 

where  A^q  is  the  water  flux  through  PEM. 

As  the  ethanol  moves  along  the  catalyst  layer,  its  flux 
decreases  due  to  ethanol  electro-oxidation  on  the  catalytic  sites. 
Thus,  the  material  balance  equation  for  ethanol  is  derived  from 
the  following  equation: 

dA^EtOH  _ _ 1 _ di 

dz  ~  12F  dz  1 


The  anode  overpotential  at  any  location  within  the  catalyst 
layer  is  defined  as: 


rj(z)  =  <ps(z)  -  <pm(z ) 


(14) 


where  <ps(z)  is  the  potential  of  the  electronic  conduction  phase  of 
the  catalyst  layer  and  (pm(z)  denotes  the  potential  of  the  ionomer 
phase.  Both  <ps(z)  and  ipm(z)  decrease  in  the  z-direction,  and  Ohm 
law  equation  for  each  phase  can  be  written  as: 


dffs 

dF 


1 

Kf 


(/  -  i) 


(15) 


d<?m  _ _ 1_. 

dz  “  Kg' 


(16) 


where  Kf  and  Kf  denote  the  effective  conductivity  of  the 
solid  phase  and  the  ionomer  phase,  respectively.  Consequently, 
the  variation  of  the  anode  overpotential  in  the  catalyst  layer  is 
given  by  the  following  equation: 


f-L  +  J ■_ 

dz  \K^  +  Kf)1  Kf 


(17) 


2.4.  Proton  exchange  membrane 


diffusion  and  hydraulic  permeation.  However,  it  was  assumed 
above  that  the  pressure  at  both  anode  and  cathode  sides  is  equal 
to  1  atm.  Therefore,  only  the  effects  of  electro-osmosis  and  dif¬ 
fusion  are  considered: 


K 1  eo  —  A^electrdrag  ~t~  A/ (jiff 


(18) 


where  Ne iectr  drag  denotes  the  water  flux  caused  by  the  electro- 
osmotic  drag,  which  at  a  constant  cell  temperature,  is  pro¬ 
portional  to  the  cell  current  density  /  and  it  can  be  expressed 
as: 

/ 

A^electrdrag  —  *M  |2()  —  (19) 

r 

where  1h2o  is  the  electro-osmotic  drag  coefficient  of  water.  /Vdiff 
is  deduced  from  the  water  concentration  gradient  through  PEM 
and  it  can  be  expressed: 


A^diff  = 


/-an  /-cath 
,vri  LH20  LH20 

un2o  7 


(20) 


where  stands  for  the  diffusion  coefficient  of  water  in  PEM, 

Zm  PEM  thickness,  CJ^q  and  C^q  the  water  concentration  at 
the  anode  and  the  cathode  side,  respectively. 

For  the  sake  of  simplicity,  based  on  the  assumption  that 
both  anode  and  cathode  are  fully  hydrated,  the  water  trans¬ 
portation  through  the  membrane  can  be  reduced  into  the  next 
equation: 


A“20  =  Ah2o- 


(21) 


Thus,  the  ethanol  flux  through  PEM  can  be  expressed  as  fol¬ 
lows: 


C' 


vEtOH 


an 
EtOH'" 


aum//tn 


c 


*vm/km  _  i 


icath 

EtOH  (,m 


(22) 


where  km  =  DgtOH//m  is  the  mass  transfer  coefficient  of  ethanol 
in  PEM,  Cg(0H  and  C|(QH  the  respective  ethanol  concentra¬ 
tion  at  the  anode  and  the  cathode  sides  of  PEM  and  i/n  = 
Mu2oNfti0/ pn^o  the  superficial  velocity  of  water  through 
PEM.  Considering  that  the  ethanol  at  the  cathode  was  from  the 
anode  due  to  crossover  and  then  it  was  oxidized  and  vapor¬ 
ized,  obviously  the  ethanol  concentration  at  the  cathode  is  much 
less  than  that  at  the  anode.  Thus,  Cj)(QH  can  be  omitted.  Conse¬ 
quently,  Eq.  (22)  is  reduced  into  Eq.  (23): 


Nm  — 
ivEtOH  — 


/'-•an  „vmjkm 
^EtOH^  m 

vm/km  v 


(23) 


This  equation  describes  the  ethanol  transportation  through  PEM, 
which  is  known  as  ethanol  crossover  and  it  will  be  further  exam¬ 
ined. 

Furthermore,  the  following  equations,  which  describe  the 
anode,  cathode  and  overall  reactions  for  a  DEFC  were  also  taken 
into  account,  while  developing  the  mathematical  model  for  the 
anode  of  a  DEFC. 


The  transportation  of  water  through  PEM  is  the  result  of 
the  combination  of  the  following  phenomena:  electro-osmosis, 


Anode  reaction  :  CH3CH2OH  +  3H20  -»  2C02  +  12H+  +  12e 

(24) 
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Cathode  reaction  :  302  +  12FI+ +  12e  ->  6H2O  (25) 

Overall  reaction  :  CH3CH2OH  +  3O2  —*■  2CO2  +  3H2O 

(26) 

3.  Results  and  discussion 

3.1.  Anode  polarization  curves  of  a  DEFC 

Fig.  2  gives  the  comparison  between  the  calculated  anode 
overpotential  using  the  present  mathematical  model  and  the 
experimental  results  [13].  The  anode  polarization  was  obtained 
by  scanning  the  potential  with  the  anode  (PtRu/C  was  the 
catalyst)  supplied  with  1.0M  ethanol  aqueous  solution  as 
the  working  electrode  and  the  cathode  fed  by  1  atm  humidi¬ 
fied  hydrogen  serving  both  as  the  counter  electrode  and  the 
dynamic  hydrogen  reference  electrode  (DHE).  The  operation 
conditions  have  been  reported  in  detail  previously  [13].  As  it 
can  be  seen  from  Fig.  2  there  is  a  good  agreement  between 
the  experimental  and  simulation  results.  This  suggests  that 
the  present  model  can  work  very  well.  It  should  be  noted 
that  the  parameters  used  in  the  present  model  were  deduced 
from  the  experimental  results  [10]  and  they  are  shown  in 
Table  1. 


Table  1 

Parameters  used  for  modeling  from  the  experimental  results  [10] 


T(°C) 

I0  (mA) 

A„  (cm2) 

Y 

30 

0.34805 

0.0285 

5.0 

0.25 

50 

0.39917 

0.0365 

5.0 

0.25 

75 

0.68331 

0.040 

5.0 

0.25 

3.2.  Effect  of  the  catalyst  layer  thickness  on  the  anode 
overpotential 

Fig.  3  shows  the  effect  of  the  catalyst  layer  thickness  on  the 
anode  overpotential.  In  Fig.  3(a),  the  effect  of  the  catalyst  layer 


Fig.  3.  Effect  of  the  catalyst  layer  thickness  on  the  anode  overpotential  at  dif¬ 
ferent  current  densities  (a)  and  different  ethanol  feed  concentrations  (b). 


662 


G.  Andreadis  et  al.  /  Journal  of  Power  Sources  157  (2006)  657-665 


Fig.  4.  Effect  of  the  catalyst  layer  thickness  on  the  ethanol  concentration  profile 
through  the  catalyst  layer,  CF,EtOH  =  2  M,  rceii  =  75  °C. 

thickness  on  the  anode  overpotential  at  different  current  den¬ 
sities  is  presented.  Two  different  ethanol  feed  concentrations 
were  considered  in  order  to  obtain  results  in  a  wide  range  of 
current  densities.  Fig.  3(b)  illustrates  the  effect  of  the  catalyst 
layer  thickness  on  the  anode  overpotential  for  different  ethanol 
feed  concentrations.  As  one  can  observe  when  the  catalyst  layer 
thickness  is  up  to  10  |xm,  the  anode  overpotential  decreases 
sharply.  This  could  be  due  to  the  insufficient  catalyst  surface 
and  thus  leading  to  the  activation  loss  [12].  Above  this  thickness 
value,  the  overpotential  continues  to  decrease,  but  slowly,  up  to 
a  certain  level  and  then  it  keeps  almost  constant  in  all  the  investi¬ 
gated  current  densities  ranging  from  50  to  1000  mA  cm-2.  More 
specifically  in  Fig.  3(a),  under  the  considered  operation  condi¬ 
tions:  /=  100mA cm-2,  CeeiOH  =  2.0 M  and  7’ceii=75°C,  the 
anode  overpotential  decreases  from  0.645  to  0.610  V  as  the  cata¬ 
lyst  layer  thickness  increases  from  10  to  20  |xm.  Moreover  for  the 
same  thickness  change,  at  /=  1000  mA  cm-2,  CF,EtOH  =  6.0M 
and  7,ceii  =  75°C,  the  overpotential  decreases  from  0.851  to 
0.849  V.  Considering  this,  the  value  of  10|xm  for  the  catalyst 
layer  thickness  was  used  in  the  following  calculations. 

The  effect  of  the  catalyst  layer  thickness  on  the  ethanol  con¬ 
centration  profile  through  the  catalyst  layer  is  given  in  Fig.  4. 
The  predicted  ethanol  concentration  through  the  catalyst  layer 
at  different  current  densities  is  presented.  It  can  be  distinguished 
that  as  the  current  density  increases,  the  predicted  ethanol  con¬ 
centration  profile  in  the  catalyst  layer  is  slightly  and  linearly 
decreased.  It  is  worthy  to  notice  that  in  the  case  of  0  mA  cm-2, 
despite  that  the  ethanol  feed  concentration  was  2.0  M,  at  the  dif¬ 
fusion  layer/catalyst  layer  interface  the  ethanol  concentration  is 
only  1 .04  M.  This  value  is  affected  mostly  by  the  mass  transfer 
coefficient  (kg)  and  less  by  the  effective  diffusion  coefficient  of 
ethanol  in  the  diffusion  layer  (77g't^j). 

eff 

3.3.  Effect  of  the  effective  protonic  conductivity  Kfl  on  the 
anode  overpotential 

Fig.  5  illustrates  the  dependence  of  the  anode  overpotential 
on  the  effective  protonic  conductivity  at  constant  operating  tern- 


Fig.  5.  Dependence  of  the  anode  overpotential  on  the  effective  protonic  conduc¬ 
tivity  at  different  current  densities,  rceii  =  75  °C.  The  inset  shows  the  dependence 
of  on  the  temperature. 

perature  (75  °C).  In  the  inset  of  Fig.  5,  the  change  of  along 
with  the  temperature  is  also  shown.  As  expected,  the  anode 
overpotential  decreases  if  the  catalyst  layer  possesses  better 
protonic  conductivity.  For  example,  at  100  mA cm-2,  when  the 
protonic  coefficient  in  the  catalyst  layer  increases  from  1  x  10-3 
to  3.2  x  10-3Scm-1,  the  corresponding  anode  overpotential 
decreases  from  0.667  to  0.648  V.  Increasing  the  protonic  coef¬ 
ficient  up  to  1  x  10_2Scm_1,  the  anode  overpotential  further 
decreases  up  to  0.64 1  V.  At  higher  current  density  values,  the  cor¬ 
responding  decrement  for  the  anode  overpotential  is  even  higher. 
Obviously,  based  on  the  above  information,  the  decreased  over¬ 
potential  can  be  achieved  by  increasing  the  protonic  conductivity 
of  the  catalyst  layer.  This  suggests  that,  in  order  to  increase  the 
effective  protonic  conductivity  in  the  catalyst  layer,  the  void 
volume  in  the  catalyst  layer  (the  space  among  catalyst  particles) 
should  be  filled  to  a  large  extent  with  ionomer  [12].  From  the 
inset  of  Fig.  5,  it  can  be  clearly  seen  that  the  effective  protonic 
conductivity  is  affected  by  the  temperature  and  it  increases  as 
the  temperature  increases.  It  is  worthy  to  note  that  the  effect  of 
the  effective  diffusion  coefficient  on  the  anode  overpotential  was 
also  investigated  and  it  was  found  that  it  had  almost  no  effect  on 
the  anode  overpotential  (Table  2). 

3.4.  Parameters  influencing  ethanol  crossover  rate 

In  Fig.  6  the  effect  of  the  current  density  on  the  ethanol 
crossover  rate  at  different  operating  temperatures  is  presented. 
As  it  can  be  seen  the  increase  of  the  current  density  affects  the 
ethanol  crossover  rate  in  two  conflicting  ways.  From  one  side, 
one  will  expect  that  increasing  the  current  density  would  def¬ 
initely  result  in  a  decrement  of  the  ethanol  crossover,  due  to 
the  fact  that  more  ethanol  will  participate  in  the  electrochemical 
reaction,  thus  decreasing  the  concentration  difference  between 
the  two  sides  of  PEM.  However,  when  high  ethanol  feed  con¬ 
centrations  (i.e.  4.0  M)  are  used,  a  volcano  behavior  is  presented 
as  shown  in  Fig.  6.  This  could  be  explained  by  the  fact  that  the 
total  ethanol  crossover  rate  is  the  combined  result  of  the  diffu- 


G.  Andre  adis  et  al.  /  Journal  of  Power  Sources  157  (2006)  657-665 


663 


Table  2 


Basic  parameter  values 


Parameter 

Value 

Ref. 

Anode  transfer  coefficient,  aa  at  75  °C 

0.04 

[13] 

Electrode  surface  (cm2) 

5 

Reference  exchange  current  density,  iQ 

0.68331 

[13] 

Order  of  reaction 

0.25 

[13] 

Catalyst  layer  thickness,  tc  (pm) 

10 

Diffusion  layer  thickness,  t&  (p,m) 

30 

PEM  thickness,  tm  (fjim) 

15 

Reference  ethanol  concentration  (M) 

0.5 

[13] 

Void  fraction  of  catalyst  layer,  ec 

0.32 

Diffusion  coefficient  of  ethanol  in  water. 

1.83  x  10“5 

[27] 

7>d,EtOH  (cm2  s  1)at75°C 

Diffusion  coefficient  of  ethanol  in  PEM, 

4.1  x  10~5 

[28] 

^EtOH  (cm”  s  1)at75°C 

Void  fraction  of  diffusion  layer,  ed 

0.8 

Protonic  conductivity  of  the  ionomer,  Km 

0.1416 

[12] 

(Scm  *) 

Effective  protonic  conductivity  in 

3.2  x  10~3 

catalyst  layer,  f  (Scm  x) 

Electronic  conductivity  of  solid  phase 

8.13  x  10-s 

[22] 

(Pt-C),  /TslScm-1) 

Effective  conductivity  of  solid  phase  in 

3.78  x  10-6 

catalyst  layer,  (Scm  1 ) 

Volume  fraction  of  ionomer  phase  in 

0.04 

catalyst  layer, 

Electro-osmotic  drag  coefficient  of 

3.16 

[12] 

water,  Xh2o 

Volume  fraction  of  solid  phase  (Pt-C)  in 

0.6 

catalyst  layer, 

Ethanol  feed  concentration  CF,EtOH  (M) 

1.0 

sion,  due  to  the  concentration  difference  between  the  anode  and 
the  cathode  and  the  electro-osmotic  drag.  As  the  current  density 
increases,  the  difference  in  ethanol  concentration  between  anode 
and  cathode  is  definitely  reduced  since  more  ethanol  is  involved 
in  the  ethanol  electro-oxidation  at  the  anode.  On  the  other  hand, 


the  electro-osmotic  drag  increases  due  to  the  fact  that  more  pro¬ 
tons  are  transported  through  the  membrane,  which  leads  to  more 
ethanol  molecules  permeated  to  the  cathode.  The  resultant  effect 
of  these  two  phenomena  can  give  a  reasonable  explanation  for 
the  observed  results.  At  lower  temperatures  (<75  °C)  such  as  30 
and  50  °C,  this  volcano  behavior  of  the  ethanol  crossover  rate 
along  with  the  current  density  is  more  obvious.  This  is  probably 
due  to  the  lower  diffusion  coefficient  of  ethanol  in  both  dif¬ 
fusion  layer  and  PEM.  Similar  predictions  have  been  obtained 
concerning  the  anode  behavior  of  a  direct  methanol  fuel  cell 
[12]. 

Fig.  7  shows  the  effect  of  the  temperature  on  the  ethanol 
crossover  rate  at  different  current  densities  and  different  ethanol 
feed  concentrations.  It  can  be  observed  that  ethanol  crossover 
rate  is  increased  with  the  temperature  increment  as  expected  due 
to  the  increased  kinetics  at  higher  temperature.  From  Fig.  7(a)  it 
can  be  seen  that  the  current  density  affects  the  ethanol  crossover 
rate  especially  at  high  current  densities.  Moreover,  the  activa¬ 
tion  energies  are  higher  at  high  current  densities.  This  could 
be  attributed  to  the  fact  that  as  the  current  density  increases, 
the  electro-osmotic  drag  is  increased  simultaneously  because 
more  protons  will  be  transported  through  PEM.  The  effect  of  the 
ethanol  feed  concentration  on  the  ethanol  crossover  rate,  at  dif¬ 
ferent  operation  temperature  is  presented  in  Fig.  7(b).  From  the 
parallel  lines,  it  is  clear  that  the  ethanol  feed  concentration  does 
not  affect  the  apparent  activation  energy  of  ethanol  crossover  in 
the  investigated  range. 

Fig.  8  illustrates  the  effect  of  the  ethanol  feed  concentration 
on  the  ethanol  crossover  rate  at  different  current  densities.  The 
ethanol  crossover  rate  increases  linearly  with  the  increment  of 
the  ethanol  feed  concentration.  In  the  case  of  CF.EtOH  =  3.0M, 
there  is  a  point  where  the  first  three  values  of  the  current  density 
are  intersected.  This  could  be  explained  by  the  volcano  behavior 
of  the  ethanol  crossover  rate  along  with  the  current  density  as 
observed  in  Fig.  6. 
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Fig.  6.  Effect  of  the  current  density  on  the  ethanol  crossover  rate  at  different  temperatures  and  different  ethanol  concentrations. 
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Fig.  7.  Effect  of  temperature  on  ethanol  crossover  rate  at  different  current  den¬ 
sities  (a)  and  different  ethanol  concentrations  (b). 

3.5.  Simulation  versus  experimental  results 

A  comparison  between  the  simulation  and  experimental 
results,  concerning  ethanol  crossover  rate,  is  illustrated  in  Fig.  9. 
The  effect  of  the  ethanol  feed  concentration  on  the  ethanol 
crossover  rate  at  different  operating  temperatures,  is  shown  in 
Fig.  9(a).  It  should  be  pointed  out  that  the  ethanol  feed  concen¬ 
tration  that  was  taken  into  account  in  the  present  results  was 


Fig.  8.  Effect  of  the  ethanol  feed  concentration  on  the  ethanol  crossover  rate  at 
different  current  densities.  Fceu  =  75  °C. 
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Fig.  9.  (a  and  b)  Comparison  of  the  effect  of  ethanol  feed  concentration  on 
ethanol  crossover  rate  at  /  =  0mAcm~2  between  the  experimental  [9,10]  and 
simulation  results. 


only  up  to  8.0 M.  There  are  more  experimental  results  concern¬ 
ing  ethanol  crossover  rate  at  higher  ethanol  feed  concentrations, 
which  will  be  reported  in  our  future  work  [10].  The  effect  of 
the  operating  temperature  on  the  ethanol  crossover  rate  along 
with  the  ethanol  concentration  is  presented  in  Fig.  9(b).  From 
the  table  shown  in  the  inset  of  Fig.  9(b),  it  is  obvious  that  there 
is  almost  the  same  apparent  activation  energy  for  all  the  investi¬ 
gated  ethanol  concentrations,  as  shown  in  Fig.  7(b).  Finally,  as  it 
can  be  seen,  there  is  a  good  agreement  between  the  experimental 
and  simulation  results,  which  further  indicates  the  validation  of 
the  present  model. 

4.  Conclusions 

In  the  present  study,  a  mathematical  model  is  developed  in 
order  to  describe  and  predict  the  behavior  of  the  anode  of  a  direct 
ethanol  proton  exchange  membrane  fuel  cell.  The  good  agree¬ 
ment  between  the  simulated  and  experimental  results  indicated 
that  the  present  developed  model  could  work  very  well.  Based 
on  the  simulation  results,  it  can  be  found  that  the  increase  of  the 
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catalyst  layer  thickness  and  the  effective  protonic  conductivity 
in  the  catalyst  layer  can  contribute  to  the  decrease  in  the  anode 
overpotential.  On  the  other  hand,  the  optimal  thickness  of  the 
catalyst  layer  is  10  p,m  considering  from  both  performance  and 
the  cost  point  of  view.  It  can  also  be  found  that  the  ethanol  feed 
concentration,  the  operation  temperature  and  the  current  density 
are  the  main  parameters  affecting  ethanol  crossover  rate  from  the 
anode  to  the  cathode  through  the  electrolyte  membrane. 
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